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Abstract. We aim to test the plausibility of a theoretical framework in which the gamma-ray 
emission detected from supernova remnants may be of hadronic origin, i.e., due to the decay 
of neutral pions produced in nuclear collisions involving relativistic nuclei. In particular, 
we investigate the effects induced by magnetic field amplification on the expected particle 
spectra, outlining a phenomenological scenario consistent with both the underlying Physics 
and the larger and larger amount of observational data provided by the present generation of 
gamma experiments, which seem to indicate rather steep spectra for the accelerated particles. 
In addition, in order to study to study how pre-supernova winds might affect the expected 
emission in this class of sources, the time-dependent gamma-ray luminosity of a remnant with 
a massive progenitor is worked out. Solid points and limitations of the proposed scenario are 
finally discussed in a critical way. 
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1 Supernova remnants, cosmic rays and gamma-ray emission 

The present generation of 7-ray telescopes, both in the GeV (Fermi, AGILE) and in the 
TeV band (HESS, VERITAS, MAGIC, CANGAROO, MILAGRO,...) is providing us with 
an unprecedented wealth of detections and observations of 7-ray-bright supernova remnants 
(SNRs). These objects have been considered for more than 70 years the main sources of 
Galactic cosmic rays (CRs) [1], but scientists are still looking for a clear-cut evidence of 
hadron acceleration in these environments. Such a smoking gun was predicted to be the 
emission of 7-rays due to the decay of neutral pions produced in nuclear collisions between 
accelerated particles and the background gas [2]. The most promising place to look at has 
been considered for some time RX J1713. 7-3946 [3, 4], the prototype of a 7-ray-bright SNR. 

Its TeV emission (between ~100 GeV and ^■lO TeV) has been extensively studied by 
several groups [5-13], who have investigated two different scenarios: besides the hadronic 
scenario outlined above, also the inverse Compton scattering of the relativistic electrons 
responsible for the X-ray synchrotron emission may produce a similar signature in the TeV 
band {leptonic scenario). Each scenario shows strong points and apparent shortcomings, 
with theorists waiting for Fermi's advent in order to discriminate between them [5-13]. The 
7-ray spectrum between ~100 MeV and ~10 GeV (hereafter the GeV band) is in fact 
expected to be significantly flatter in the leptonic scenario with respect to the hadronic one. 
Very recently Fermi-LAT eventually showed with unprecedented resolution the morphology 
of the 7-ray emission in the region around RX J1713. 7-3946, favouring a leptonic rather than 
a hadronic emission from the SNR shell and eventually undermining physicists' hope to have 
found the smoking-gun for hadron acceleration in SNRs [4]. 

We start our analysis by building up the list of Galactic SNRs detected in 7-rays, 
followed by the list of other possible candidates (table 1). This latter class of objects is 
made of sources whose association with a SNR is plausible but not fully confirmed because 
of the tangled morphology of the emitting region and/or because of the "pollution" due to 
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other potential 7-ray sources, like pulsars and pulsar wind nebulae. Generally speaking, a 
clear detection of a SNR shell is an exceptional fact (probably the best example is Vela 
Jr, with SN1006 also showing a sharp bilateral emission consistent with its non-thermal X- 
ray morphology). In addition, it is necessary to remember that in many cases the angular 
resolution of the telescopes does not allow for a detailed morphological study of extended 
sources. The association of the 7-ray emission with molecular clouds (MCs) is also very 
common: in table 2 of ref. [14] the reader can find a thorough collection of Galactic SNRs 
known to be in physical contact with MCs, i.e. showing evidences of interaction between the 
forward shock and the MC. Very generally, also the 7-ray emission detected from sources 
not listed in [14] may be (and has often been interpreted to be) due to MCs surrounding the 
SNR rather than to the shell itself. 

In the present work we consider only particles accelerated at SNR forward shocks, 
neglecting both the difficulties related to particle acceleration in dense, partially neutral, 
environments and to the possibility for the emission to come from MCs illuminated by CRs 
escaping the SNR. We will comment more on these points below, in section 5. This approach is 
motivated by the fact that basically all of the 7-ray-bright SNRs show quite similar properties, 
and in particular energy spectra steeper than E~'^, independently from their association with 
MCs. Within a hadronic scenario for the 7-ray emission, the photon spectrum above ~100 
MeV has to be parallel to the proton one, therefore the observational evidences reported 
in table 1 can put a strict constraint on the physical mechanisms responsible for hadron 
acceleration in this class of sources. 

The success of the so-called SNR paradigm for the origin of Galactic CRs, besides 
the energetic argument already put forward by Baade and Zwicky [1], mostly relies on the 
generality of the spectrum of accelerated particles produced via Fermi's first order mechanism 
at strong shocks. In fact, during the late '70s several authors independently realized that 
diffusive shock acceleration (DSA) leads to power-law spectra 

iV(S)ocE-"; a = ^^, (1.1) 

r — 1 

whose slope does not depend on the microphysical details of the scattering process but only 
on the shock compression ratio r, i.e. the ratio between the density of the hotter, shocked 
plasma and of the colder, unperturbed one [15-20]. In particular, for strong shocks (namely 
with sonic Mach number much larger than 1) we have 

r = (1.2) 

7-1 

and hence, for standard monoatomic gas with adiabatic index 7 = 5/3, we get N{E) oc E~'^. 
Such a power-law spectrum is particularly appealing since it is in decent agreement both 
with the multi-wavelentgh observations of SNRs and, once corrected for propagation in the 
Milky Way, with the diffuse spectrum of Galactic CRs measured at Earth as well. 

Nevertheless, as soon as the first quantitative calculations about DSA efficiency were 
carried out, people realized that such a process may channel a large fraction (even more than 
90%) of the fluid ram pressure into CRs. In such a scenario, accelerated particles can no 
longer be viewed as test-particles and the CR population has to be treated as an additional 
component entering the hydrodynamical equations for conservation of mass, momentum and 
energy [21, 22]. This two-fluids approach to the non-linear theory of DSA (NLDSA) has soon 
been followed by kinetic approaches to the problem, in which also the information about the 
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momentum distribution of accelerated particles is retained. In this case the equations of the 
shock hydrodynamics, modified by the inclusion of CRs, are solved along with a description of 
the CR transport, typically accounted for via a diffusion-convection equation or a numerical 
(Monte Carlo) treatment. Comprehensive reviews on CR- modified shocks can be found 
in refs. [23-26], while in ref. [27] different kinetic approaches to NLDSA (semi-analytical, 
numerical and Monte Carlo) are summarized and compared. 

The most important prediction of any NLDSA theory is that the pressure in accelerated 
particles around the shock leads to the formation of a precursor in which the upstream 
fluid is slowed down and compressed. As a consequence, the proper shock becomes weaker 
(it is in fact referred to as suh shock) and provides a reduced heating of the downstream 
plasma. Accelerated particles diffusing in such a modified fluid profile feel different (averaged) 
compression ratios and, since high-energy particles have diffusion lengths larger than low- 
energy ones, the overall spectrum can no longer be a simple power-law. High-energy particles 
probe in fact the whole precursor, and hence feel the total compression ratio -Rtot > 4, while 
low-energy particles are confined close the subshock, whose compression ratio is instead 
-Rsub < 4. As a consequence of the intrinsic property of Fermi's first order mechanism — see 
eq. (1.1) — the resulting spectrum is predicted to be concave, i.e. steeper (flatter) than 
at low (high) energies. 

This effect for high-energy particles (above some GeV, as a rule of thumb) can be also 
accounted for as a pure hydro dynamical effect: the contribution to pressure and energy by 
relativistic particles (whose adiabatic index is 7 = 4/3) actually makes the total (gas -|- CRs) 
fluid more compressible. The global effective adiabatic index can be written (see e.g. [28, 29]) 
as: 

15 + 36^ 

7cff = :t -1 . ^ , (1-3) 

where ^cr is the fraction of the bulk pressure converted in accelerated particles at the shock. 
Since 4/3 < 7cfr < 5/3, we have 4 < r^^ ^ 7 and finally 1.5 < a < 2: the more efficient 
the acceleration, the flatter the spectrum of the accelerated particles. Such a simple trend is 
recovered in all of the NLDSA models and becomes more and more marked if also particle 
escape is accounted for: in this case, in fact, the shock behaves as partially radiative and the 
compression ratio might become much larger than 7, eventually leading to spectra as flat as 
~ at the highest energies (see e.g. [26, 30]). 

Since the spectrum of hadronic 7-rays maps the spectrum of accelerated particles un- 
equivocally, concave photon spectra would naturally represent the smoking gun for very ef- 
ficient hadron acceleration in SNRs. In addition, NLDSA naturally predicts pretty hard 
spectra at large energies, so that hadronic 7-ray emission in the TeV range should typically 
show a spectrum harder than E~'^ . Neither the former, stronger evidence, nor the latter, 
more general, one finds any support in the observations summarized in table 1: does this fact 
imply that the 7-ray emission can not be of hadronic origin? Or that SNRs are not efficient 
CR accelerators? Or, even, both of them? 

In the following sections we show how including an additional piece of information, 
namely the amplification of magnetic field due to plasma instabilities excited by the streaming 
of relativistic particles, it is possible to outline a physically consistent scenario in which SNRs 
are efficient factories of Galactic CRs and the observed 7-ray emission may nevertheless be 
explained within a hadronic scenario. 

In section 2 we introduce a phenomenological model accounting for particle acceleration 
and magnetic field amplification at shocks. In section 3 we outline a quite general treatment 
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of the SNR hydro dynamical evolution, also allowing for the presence of pre-SN winds con- 
tributing to create a non-trivial circumstellar environment the shock propagates into. We 
then consider in section 4 the case of remnant with a massive progenitor, in order to illustrate 
the effects induced by the winds launched during pre-SN stages on the time evolution of the 
7-ray emission, to some extent generalizing the pioneering work in ref. [2]. In particular, 
such a choice allows us non only to assess the fundamental role of magnetic field amplifi- 
cation in producing rather steep spectra of accelerated particles, but also to argue that the 
age distribution of 7-ray-bright SNRs correlates better with core-collapse SNRs rather than 
with type la SNRs expanding into homogeneous environments (section 4.1). Our findings 
are finally commented in section 5, where we critically discuss the modeling of magnetic field 
amplification and the possible role of MCs often surrounding SNRs. We conclude in section 
6. 

2 Modeling the spectrum of cosmic rays 

During the last decade. X-ray observations of young SNRs have shown evidences of bright 
narrow rims of non-thermal origin, which have been interpreted as due to synchrotron emis- 
sion from ultra-relativistic electrons. Moreover, from the thickness of these rims it has been 
possible to put a lower limit for the magnetic field immediately behind the shock, finding 
evidences of magnetic fields as high as a few hundred ^G, almost two order of magnitudes 
larger than the standard interstellar field (see e.g. refs. [62, 82-85]). 

These amplified magnetic fields are very likely produced by the super- Alfvenic streaming 
of particles accelerated at the shock via plasma instabilities [19, 86-95]. In such a scenario, 
CRs are able to generate by themselves the magnetic turbulence responsible for their own 
diffusion in a non-linear interplay which eventually allows accelerated particles to reach en- 
ergies consistent with the steepening (the knee, around 3 x lO^GeV) observed in the diffuse 
spectrum of Galactic CRs and to explain it as due to the maximum energy achievable in 
SNRs [96]. An evidence that magnetic field amplification occurs upstream comes from the 
narrow extension of the X-ray emitting region ahead of the shock in SN1006, as it has been 
put forward by exploiting very detailed Chandra maps [97, 98]. 

Two main categories of streaming instabilities are usually accounted for: resonant and 
non-resonant ones, according to the relation between the Larmor radius of the relativistic 
particles and the wavelength of the excited modes. Resonant streaming instability has been 
known since the '70s [19, 86-88], while more recently T. Bell worked out a class of short- 
wavelength modes whose growth may be even faster [89, 90]. A kinetic approach to streaming 
instability can show how resonant and non-resonant modes are two faces of the same coin 
and that non-resonant modes grow faster than resonant ones only for large shock velocities 
[91]. For this reason, and also because the role of non-resonant modes in CR diffusion is 
still to be addressed, we focus our attention only on resonant modes, checking a posteriori 
the possible role of Bell's modes. In addition, we bypass the time-dependent details of the 
magnetic turbulence growth, assuming that field amplification proceeds until saturation is 
achieved. This should happen when the normalized pressure in the amplified magnetic field 
is (see [30, 99] and references therein) 

^ = « (2 1) 

PVI 8npVl 2Ma' ^ ■ ^ 

with Ma = Vsh/vA_ the Alfvenic Mach number. It is important to stress that the saturation 
above is calculated as a non-linear extrapolation of the result obtained within a quasi-linear 
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theory for magnetic field amplification, where Ma is taken in the background field, Bq. 
However, since a configuration with perturbations 5B 3> Bq is rather peculiar and likely 
unstable (also because 5B _L Bq), we could phenomenologically assume that the Alfvenic 
Mach number entering eq. (2.1) should be taken as something closer to Ma{B) rather than 
to Ma{Bq). This choice would also imply that the Alfven velocity va = B/ yJAirp associated 
with the turbulence should be calculated in the amplified magnetic field as well. Finally, 
since CRs actually scatter against magnetic irregularities rather than against the fluid itself, 
when VA is not negligible with respect to fluid velocity u, the compression ratio actually felt 
by accelerated particles has to be accordingly modified, namely: 



U2 



Ul + VA,l 
U2 + VA,2 ' 



(2.2) 



where subscripts 1 and 2 refer to pre- and post-shock quantities, respectively. Since Alfven 
waves are generated upstream by streaming instability, they must travel against the fluid (in 
direction opposite to the CR pressure gradient), while it seems reasonable to assume isotropic 
turbulence downstream, which implies va,2 = 0, on average. Moreover, when CR acceleration 
is efficient Ui = ui/Vsh — 1 — ^cr and the equation for streaming instability saturation reads 
(see eq. 42 of [99]) 



Pb,i 



-3/2 



1 



from which 



Ma,: 



2 {l-icr) 



AMa, 



5/2 



(2.3) 



(2.4) 



Ccr 2 

Finally, dividing numerator and denominator in eq. (2.2) by ui and retaining the general 
formula for the shock compression ratio we get 



1 



1 



Ma, 



7efr + 1 



7eff - 1 + 2/M2 



1 



Ccr(2 ^cr) 
2(1 -Ccr.)^/^ 



(2.5) 



where Ms = V^h/cs is the sonic Mach number of the shock and 7cfr is defined in eq. (1.3). 

In this work we consider the CR acceleration efficiency ^cr as given, keeping in mind 
that, in order to account for the energetics in Galactic CRs, a fraction of about 10% of the 
SN kinetic energy has to be channelled into accelerated particles (see e.g. the review in 
ref. [100]). The spectral slope of the accelerated particles is thus simply given, as in eq. (1.1), 

by 

(2^6) 



a 



1 



In figure 1 we show the calculated spectral index a (dark lines, left axis) and the related 
amplified magnetic field upstream of the shock, Bi (light lines, right function of the 

CR acceleration efficiency ^cr, taken to vary between zero and 20%. We fixed temperature 
and density of the background gas in Tq = 10^ K and po = O.lm^cm"^, hence the dependence 
of a on Ms is accounted for showing different curves with different shock velocities (as in the 
legend), corresponding to Ms ~ 50,17 and 5, respectively. 

From figure 1 it is clear that for large Mach numbers a non-negligible CR efficiency is 
required to produce particle spectra steeper than E~'^ while, on the other hand, for sufficiently 
low Ms the standard reduction of the compression ratio predicted for non-strong shocks is 
amplified even further with increasing ^cr- Moreover, we notice that the effect of taking 
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Figure 1. Slope of the energy spectrum of accelerated particles (dark lines, laft axis) and upstream 
amplified magnetic field (light lines, right axis) as a function of the CR acceleration efficiency (see 
eqs. 2.3-2.6). Temperature and density of the background plasma are fixed in Tq = 10^ K and 
Pa = 0.l7TipCm^'^, while different lines correspond to different shock velocities, as in the legend. 



the saturation of the instabihty in the amplified field rather than in the background one 
completely suppresses the spectral flattening usually predicted by NLDSA theory, so that r 
never grows beyond 4, neither for large ^cr- This prediction has however to be taken with 
a grain of salt, since the simple model outlined here does not account for any concavity in 
the spectrum of accelerated particles and therefore can be viewed as a heuristic approach 
which loses its validity whenever a very marked precursors is produced by a very efficient CR 
acceleration. However, since at the moment there are no observational evidences that very 
large compression ratios are achieved in SNR environments (the most clear example is likely 
represented by Tycho's SNR, where 4 < r < 7 is inferred [101]), there are neither consistency 
issues or observational shortcomings with our working hypotheses. 

Finally, we also notice that the amplified magnetic field so produced (right axis in 
figure 1) are in decent agreement with the ones needed to account for the measured thickness 
of non-thermal X-ray rims in several young remnants, which require typical downstream 
magnetic fields B2 — rBi ~ 50 — 300/iG [83, 84, 102]. Moreover, when instability saturation 
is given by eq. 2.3, any dependence on the strength of background magnetic field Bq is washed 
away: the only necessary requirement is the presence of a non-vanshing component of the 
background field parallel to the CR pressure gradient, i.e. parallel to the shock normal. 

It is worth stressing that the possible role of a finite scattering center velocity has 
been put forward since when DSA model were first proposed [19], but only in very recent 
times, motivated by the discovery of efficient magnetic field amplification in young SNRs, 
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this additional piece of information has been accounted for in order to explain the observed 
phenomenology (see e.g. [99, 103-105]). Here we show, for the first time as far as we know, 
the crucial implication of this phenomenon in the calculation of the hadronic emission from 
SNRs: the inclusion of such an effect, in fact, at the moment seems to be the only reasonable 
way to obtain efficient particle acceleration along with steep spectra. 

3 Modeling the target: non-homogeneous circumstellar environments 

Modeling particle acceleration is only half of the problem of understanding the properties of 
7-ray-bright SNRs, the remaining half being the environment they expand into. This latter 
ingredient is important because regulates the number density of accelerated particles and 
provides targets for p-p scattering, as well. The radiative signature in 7-rays due to shock 
propagation into the dense (n ~ 0.1 — lcm~^) homogeneous ISM has been already extensively 
investigated in ref. [2] and in papers based on it. Here we want to set up a more general 
scheme able to account also for complex circumstellar environments, like the ones in which 
SNRs with massive progenitors expand into. 

Core-collapse SNe are the final evolutionary step of very massive stars, which enter 
the main sequence (Hydrogen burning stage) with more than 5-10 solar masses (Mq) and 
progressively lose large fractions of their mass in the shape of stellar winds. In general, 
the calculation of the exact amount of mass which goes into these winds, and hence of the 
remainder which goes into the SN ejecta or into the central compact object, depends on many 
variables, the most important of which are likely the initial mass and the metallicity. Let us 
consider now two major kinds of stellar mass ejection which may occur during pre-SN stages: 
red-giant and Wolf-Rayet winds. 

A single massive star may produce one or even both of them during its life, thus we 
consider here a quite general case in which a cold, slow and dense red-giant wind is launched, 
followed by a hot and fast Wolf-Rayet one. The latter is expected to penetrate through the 
former, eventually excavating a large, rarefied, cavity around the star, often referred to as a 
hot bubble. The fine structure of such an interplay is outlined, e.g., in refs. [106-108], but for 
our purposes it is worth sketching a rather simpler picture, similar to the one illustrated in 
figure 3 of ref. [109]. 

The circumstellar medium is hence taken as spherically symmetrical and as in the fol- 
lowing. The innermost region is occupied by the dense and cold red-giant wind, whose density 
profile depends on the radial coordinate r as 

= Tir^ = 3.5mp— ^ — cm-3 , 3.1 

where rup is the proton mass, M = M_5,0lO-^M0/yr is the mass-loss rate during the wind 
activity, Vw = T4,6l0^cm/s is the wind velocity; both quantities are normalized to commonly 
accepted values. The wind temperature is expected to be T = 10^ — 10^ K [109]. The 
magnetic field in such a wind is usually thought to be in the shape of a Parker spiral, 
and therefore exactly perpendicular to the shock normal: this fact may be relevant for our 
purposes since without a magnetic field component parallel to the CR gradient, the standard 
streaming instability introduced above should not be very effective. On the other hand, it 
is very likely for the hydro instabilities due to the Wolf-Rayet wind crossing to produce a 
turbulent (disordered) magnetic field: as a consequence the shock should be predominantly 
oblique rather than perpendicular and therefore all of the arguments about magnetic field 
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amplification put forward in section 2 should qualitatively apply to this region as well. The 
radial extension of the wind zone Ry^ is given by the total mass ejected during the red- giant 
stage, Mw, i.e.: 

fRw ]\/f T/ „ 

= / 47rp^(r)dr = 1.4 ^"'"'® '^"'^ pc . (3.2) 

Jo M_5^Q 

Beyond R^ we find the hot bubble excavated by the Wolf-Rayet wind, which is usually taken 
as much more homogeneous and rarefied. Fiducial values for temperature and density are 
Ti) ~ 10^ K and pb = O.OlmpCm"^. Like for the dense wind, the extension of the hot bubble 
can be determined by the total mass ejected during the Wolf-Rayet stage, namely: 

Mb = J 47rpb(r)dr ^ fib = 10.23 |^Afb,0 ' ^ " j pc. (3.3) 

Phenomena occurring at the transition between the hot bubble and the interstellar medium 
(ISM) are thoroughly studied in ref. [109] but, for our purposes, we can simply assume a 
sharp discontinuity by posing, for r > R^, Tq = 10*^ K, po/nip = lcm~^. 

Massive stars typically lose a large fraction of their initial mass in stellar winds, there- 
fore only a minor fraction goes into the ejecta during the SN explosion. In terms of SNR 
evolution, this fact implies that the remnant of a core-collapse SN enters its Sedov- Taylor 
stage when its forward shock is still propagating inside the progenitor wind, and precisely 
at RsT — l-4Mej^0pc for the case above. For smaller radii (i.e. during the ejecta-dominated 
stage), position and velocity of the forward shock can be calculated with suitable self-similar 
analytical solutions [110]. Here we adopt the same recipe as in eq. 15 of ref. [Ill], which 
reads: 

^ 99R%1 (-^^ 1 yr (3.4) 




Kh(/?sh) ^ 8800R:^](J^ ( -f^^^ ] km s-\ (3.5) 



For r > RsT the forward shock propagation into the complex circumstellar environment 
can be, instead, calculated by adopting the so-called thin-shell approximation, i.e. neglecting 
the spatial structure inside the shock and deriving position (i?sh) and velocity (Fsh) of the 
forward shock by assuming the total mass (ejecta -|- swept-up) to be concentrated in a thin 
shell around -Rgh- The shock evolution is hence obtained by solving the time-dependent 
equations for continuity of mass, momentum and energy and reads (see appendix of ref. [Ill] 
and ref. [112]): 

rRsh 

M(i?,h) = Mej + 47r / drr^r); (3.6) 
Jo 

£:(i?,h) = £sN - 47r / * drr^T^scir); (3.7) 
Jo 

i(fish) = TF^; A = (3.8) 







Kh(-Rsh) 



7cff + 1 



2A /■^* 



M^R,^)R^ 



— / drr^~^E{r)M{r] 
sh -^0 



1/2 



(3.9) 
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where Mej and SsTSS = "fsilO erg are respectively the mass and the kinetic energy of the SN 
ejecta. Tesc is the flux of energy carried away from the system by escaping CRs, assumed to be 
negligible during the ejecta-dominated stage (see ref. [113] for details). This non-adiabatic 
term adds up to the standard approach of including the CR contribution by adopting a 
suitable adiabatic index for the fluid, as in eq. (1.3), and allows us to test an additional effect 
of an efficient CR production on the shock dynamics (also see ref. [114]). The evolution of 
the remnant is followed until the end of the Sedov- Taylor stage. A generalization of this 
method able to account for the transition to the radiative stage could be performed following 
the analytic approach put forward in ref. [115], but it is beyond the main goal of this paper. 

4 What observations say 

The mathematical apparatus put forward in the previous sections may be applied to basically 
any given SNR since it can easily account either for red-giant or Wolf-Rayet wind only, or 
even describe the evolution of type la SNRs expanding in the dense and cold homogeneous 
ISM. The CR acceleration and the argument about the role of the amplified magnetic field 
hold basically independently of the circumstellar environment. 

In this section we choose to focus our attention on remnants produced by core-collapse 
SNe for several reasons: 1) core-collapse (type Ib/c and type II) SNe are more common than 
type la ones, accounting for about 70-80% of total Galactic SNe [116]; they are therefore 
more representative for our statistical analysis; 2) in young SNRs (with strong shocks) the 
spectral slope depends only on the CR acceleration efficiency, as showed in figure 1, therefore 
the key point concerning the need to explain the steep spectra observed is almost independent 
of the details of the circumstellar medium; 3) it is indeed interesting to investigate the time 
evolution of the 7-ray luminosity when the shock propagates into a circumstellar medium 
significantly modified by the stellar winds launched during pre-SN stages. This latter aspect 
is complementary to the analysis already led for a homogeneous circumstellar medium ([2] and 
papers based on it) and may give an important insight into the nature of the progenitors of 7- 
ray-bright SNRs. In particular, in section 4.1 we point out how the observed age distribution 
of the sources in table 1 is more consistent with the one expected from SNRs whose progenitor 
blew a rarefied hot bubble rather than from SNRs expanding in the unperturbed ISM. 

Let us consider a massive star whose total mass in ejecta, red-giant wind and Wolf-Rayet 
wind (i.e. the hot bubble mass) are respectively Mej = 4:Mq, = 12Mq and = 16Mq. 
Let us also take standard values £51 = V^fi = M-^^q = 1. In addition, CR acceleration 
efficiency is fixed in = 0.1 and the instantaneous flux of energy carried away by escaping 
CRs Jcsc is self-consistently calculated according to eq. 24 in ref. [113]. 

The evolution of relevant physical quantities is depicted in figure 2. The vertical solid 
line marks the boundary between the ejecta-dominated and the Sedov- Taylor stages, while 
vertical dashed lines indicate the SNR evolutionary stages in terms of the medium the forward 
shock propagate into: from inside to outside we have the red-giant wind, the hot bubble and 
finally the ordinary ISM, as described above. 

It is interesting to notice that the standard adiabatic solution for the shock position and 
velocity is recovered until the amount of energy carried away by escaping particles (around 
10% of £sn) begins to play a non-neglibile role in the shock dynamics, as showed by the 
steepening of the velocity curve after ~ 3000 yr. The SNR radius eventually almost stands 
when the forward shock encounters the dense and cold ISM, slowing down quite abruptly. 
Another interesting evolution worth noticing is the one of the downstream magnetic field 
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Figure 2. Time evolution of relevant physical quantities for a fixed CR acceleration efficiency ^cr = 
0.1. The vertical solid line indicates the transition between ejecta-doniinated and Sedov- Taylor stages, 
while vertical dashed lines, from left to right, mark the boundaries of wind zone, hot bubble and ISM, 
as in the labels (see also the description in section 3). 

{B2), which is efficiently amplified by CR streaming instability during early stages and then 
drops to its typical interstellar value (a few ^G) as the shock velocity decreases. Such a 
trend is in qualitatively good agreement with the observational evidence that only young 
SNRs are very bright in non-thermal X-rays, and it is also in quantitative agreement with 
the strength — a few hundreds /iG — of the downstream magnetic fields inferred in young 
shell-type SNRs [62, 82-85]. 

In figure 3 the predictions for the slope of the CR spectrum, a, (and therefore for the 
photon index above ~100 MeV in a hadronic scenario) are plotted as a function of the SNR 
age for three different acceleration efficiencies: ^cr = 0.01,0.1 and 0.2, as in the legend. The 
upper panel refers to observation by 7-ray satellites in the 100 MeV-50 GeV energy range, 
while the bottom panel illustrates data from Cherenkov telescopes in the region between 50 
GeV and 50 TeV. Source data collected in table 1 are divided into three categories: SNRs 
(black squares), SNR candidates (red diamonds) and finally sources showing evidences for a 
cut-off around 1-10 GeV (green triangles). 

The most striking result of the present calculation is that magnetic field amplification 
due to CR streaming instability may help to explain the spectral slopes inferred by 7-ray ob- 
servation of SNRs, provided the acceleration efficiency to be as large as 10-20%. In fact most 
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Figure 3. 7-ray-bright SNRs detected in the GeV (top panel) and in the TeV (bottom panel) band 
as in table 1. Different lines correspond to the CR spectral slope as a function of time, for different 
acceleration efficiencies as in the legend. 
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of the data points, both in the GeV and in the TeV range, fah above the hne corresponding 
to the inefficient case (^cr = 0.01, dashed hne) and he between the moderately efficient case 
(?cr = 0.1, dot-dashed hne) and the efficient one (^cr = 0.2, solid hne). More precisely, the 
dependence of the spectral slope on the CR acceleration efficiency is rather strong above 
^cr ~ 0.1, so that efficiencies larger than 20% are never required. This fact is particularly 
appealing also because NLDSA has widely shown how intrinsically difficult it is to produce 
spectra much steeper than E~'^ and, at the same time, to achieve large acceleration effi- 
ciencies. It is important to stress that results here have been worked out not within a fully 
non-linear theory of DSA, but rather within a hybrid approach which accounts in a decent 
way for the CR feedback on the SNR evolution. Such a treatment neither retains any kinetic 
information about the spectral concavity, nor evaluates in a self-consistent way the actual 
position of the cut-off in the proton spectrum. In this respect, it is indeed useful to recall 
that both the preliminary NLDSA calculations in section 5.1 of ref. [99] and the more com- 
plete ones outlined in ref. [117] demonstrate that it is possible to get a moderately efficient 
scenario consistent with CR spectra significantly steeper than the test-particle prediction by 
accounting for an effective scattering center velocity, as in section 2. As a general considera- 
tion, as long as we are dealing with efficiencies of order ~ 10%, non-linear effects like a strong 
upstream precursor implying a strong spectral concavity can be safely neglected (e.g. [105]), 
in turn justifying a posteriori the assumptions made in the present work (also see section 5.1 
for further comments on this point). 

The simple approach proposed here includes many physical ingredients which allow us 
not only to study in a consistent way the remnant evolution, but also to have an estimate of 
the properties of the expected 7-ray-emission. A prediction of this model is that, in young 
SNRs, the CR spectral index may depend only on the CR acceleration efficiency. In fact, as 
long as both sonic and Alfvenic Mach numbers are large, a depends only on Ma,i which is, in 
turn, a function of ^cr only (eq. 2.4). It is worth stressing that this fact is a direct consequence 
of having assumed the Alfven velocity relevant for streaming instability saturation to be the 
one in the amplified magnetic field, and not in the background one. 

When Ma decreases because of the low density (as in the hot bubble) or when the 
shock slows down because of the swept-up mass, the spectral slope gradually increases: this 
is the reason why middle-age and old SNRs are expected to show rather steeper spectra. In 
addition, also the energy carried away by escaping CRs may be relevant for the determination 
of the spectral slope: lines corresponding to different ^cr in figure 3 are in fact parallel to 
each other until f^esc becomes a non-negligible fraction of £sn- This energy depletion due to 
CR escape might lead, in principle, to arbitrarily steep spectra for old SNRs, as suggested 
by solid lines in figure 3, and/or to an early death of the remnant in terms of its non-thermal 
activity. 

Some TeV data points seem to fall quite outside the region consistent with our pre- 
dictions, but most of these correspond to sources showing evidence for a cut-off in the 7- 
ray-emission around 1-10 GeV (green triangles in figure 3). For these sources the physical 
scenario depicted above does not apply because their TeV spectrum may resemble more an 
exponential tail rather than an ordinary power-law. Such a warning also applies, to minor 
extent, to basically all the TeV data points since, apart from few well-studied SNRs show- 
ing clear-cut evidences for a cut-off around some TeV, the superposition of multiple sources 
and/or the low statistics do not allow us to distinguish between a cut-off spectrum or a steeper 
power-law. In this respect, the wealth of new, high-resolution data obtained with present and 
future 7-ray experiments is expected to provide more and more accurate information about 
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Figure 4. Time evolution of SNR 7-ray luminosity, both in the GeV (upper lines) and in the TeV band 
(lower lines). Solid, dashed and dot-dashed lines correspond to different CR acceleration efficiency, 
as in the legend. Vertical lines illustrate the evolutionary stages for the case ^cr = 0.1, as in figure 2. 

the statistics of SNR 7-ray spectra, in turn allowing a better and better comprehension of 
the relationship among remnant age and both spectral slope and instantaneous maximum 
energy Emax of accelerated hadrons. 

An interesting information which could be extracted from 7-ray spectra, when inter- 
preted to be of hadronic origin, could in fact be the distribution of Emax as a function of 
the SNR evolutionary stage. While there is a wide consensus on the expected increase of 
the age-dominated Emax during the ejecta-dominated stage (see e.g. refs. [96, 118, 119]), 
the actual behavior of Emax{t) during later stages is still an open question. The expected 
decrease of the magnetic turbulence level suggests a space- limited Emax , determined by a less 
and less efficient confinement of high-energy CRs. Unfortunately, quantitative calculations 
of the actual magnetic field damping in middle-age SNRs depend on a plethora of plasma 
processes very difficult to keep under control from ffist principles [111, 120]. 

4.1 Evolution of the gamma-ray emission 

In this section we work out the time evolution of the expected 7-ray luminosity by adopting 
the same techniques of ref. [2] , not in the case of the homogeneous ISM but for the circum- 
stellar profile outlined above, instead. At any given time, in fact, it is possible to estimate 
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the total number of expected 7-rays with energy above Eq as 



W^{> Eo,t):^q^{Eo,a)M{t) 



(4.1) 




where q^{EQ,a) contains the whole information about nuclear interactions and about the 
shape of the proton spectrum (see table 1 in ref. [2]), M{t) is, as usual, the total swept-up 
mass and Scrit) is the total energy in CRs, i.e.: 



In figure 4 the total luminosity in photons per second above 100 MeV (upper lines) and 
above 1 TeV (lower lines) is plotted; different lines correspond to the three CR acceleration 
efficiencies considered above. The most interesting feature is that there is a minimum of 
the emission when the forward shock propagates into the hot bubble, which can easily be 
interpreted as a consequence of the scarce amount of targets available during this stage. The 
total luminosity may be orders of magnitudes larger at earlier and later stages, inside the 
dense wind and inside the ISM, respectively. In passing by, we also notice that in the efficient 
case (^cr = 0.2, solid line in figure 4) the impact of the forward shock with the dense ISM 
occurs at later times with respect to less efficient cases: this fact is a consequence of the 
energy carried away by escaping CRs, which efficiently slows down the shock already during 
its propagation inside the hot bubble (compare the position of steep rise of the solid curve 
around 10^ yr with respect to the vertical dashed line marking the transition between the 
bubble and the ISM for the case ^cr = 0.1, as in figure 2). 

By looking at the inferred ages of 7-ray-bright SNRs, it is possible to notice a sort of 
bimodal distribution: most of the SNRs seem to be either quite young (less than 3000 yr-old) 
or rather old (more than 10000 yr-old). The statistics may still be quite low to claim a solid 
evidence for a dichotomy in the distribution but, on the other hand, we showed that it is 
quite natural to expect a significant reduction of the hadronic 7-ray-emission when the shock 
is propagating inside the hot bubble. For massive progenitors, this SNR stage should begin 
around a few thousands yr after the SN explosion and last as much time, depending on the 
mass-loss rate and total mass ejected during the red-giant and the Wolf-Rayet wind periods. 

Since in the Milky Way core-collapse SNe are more frequent by a factor 4-5 [116] than 
type la ones, which explode in the homogeneous ISM and are not expected to show a drop 
in the 7-ray luminosity at intermediate stages, it is plausible for the signature of a bimodal 
distribution to pop up in a representative sample of Galactic 7-ray-bright SNRs. 

On the other hand, SNRs with type la progenitors show a peak in the emission around 
5000 yr [2], exactly in the less populated region of figure 3. Therefore, the age distribution of 
7-ray-bright SNRs listed in table 1 is quite at odds with a scenario accounting only for SNRs 
expanding into the unperturbed ISM and seems to favor the correlation with SNR expanding, 
for a limited stage of their life, in the rarefied bubble excavated by fast pre-SN winds. In 
addition, the presence of dense red-giant winds may help to explain the pretty large number 
of very young SNRs (say with ages below 3000 yr) detected in 7-rays, which would not be 
accounted for if young SNRs expanded in the hot bubble only. The amounts of mass going 
into the different winds in the present example (figure 2) have hence been chosen in order 
to match the stage of propagation in the bubble in correspondence with the underpopulated 
age region in figure 3 (between 3000 and 10000 yr). 




(4.2) 
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It is however important to stress that core-cohapse SNe are typicahy found in cluster 
of young stars and are thus often associated with MCs which may significantly enhance the 
7-ray luminosity by providing extra targets. If this were the most common case among 
the objects in table 3, however, it would not be easy to unequivocally associate the age 
distribution of 7-ray-bright SNRs with the nature of their progenitor. 

This topic is indeed particularly intriguing since the actual factories of Galactic CRs 
have been preferentially related either to isolated SNRs exploding in the standard ISM on 
the basis of the CR chemical composition (see e.g. refs. [121]), or to SN-rich superbubbles 
[122]. Future 7-ray observations might unravel the question, providing an interesting insight 
into the nature of the medium CR accelerators explode into by identifying possible "gaps" 
in the number of detected sources as a function of their age. 

4.2 GeV/TeV connection 

As long as a SNR is observed in a single band, it is very difficult to unambiguously identify the 
mechanisms responsible for its non-thermal emission and, in fact, every studied SNR shows a 
7-ray emission that could be (and has been) fitted either within a leptonic or within a hadronic 
scenario. The idea behind this point is that a spectrum steeper than E~'^ can be accounted for 
either as the result of the decay of neutral pions produced in hadronic interactions or as due 
to the relativistic bremsstrahlung from accelerated electrons. In addition, in many cases also 
the inverse-Compton scattering of relativistic electrons on a suitable photon background can 
provide a relevant contribution in the same energy region. Inverse-Compton photon spectra 
produced by electrons distributed as a power-law E~'^ are still power-laws, with spectral 
index f3 = {a — l)/2 < a. However, when observed in a limited range of energies close to 
or above the cut-off, also a rather flat inverse-Compton contribution may mimic a steeper 
spectrum, hence it is very hard to rule out such an emission mechanism even if basically all 
of the observed spectra are steeper than (see e.g. [123] for a wider discussion). 

One could expect such an ambiguity to be washed away when broadband (i.e. from 
sub-GeV to multi TeV) data were available, because of the information enclosed in the 
relative normalization of the GeV and TeV spectra. This makes perfect sense, but at the 
moment there are only a few SNRs detected in both bands (see table 1): in Cassiopeia A the 
superposition of relativistic bremsstrahlung and inverse-Compton contributions provides a 
7-ray spectrum consistent with the data as well as the hadronic scenario [45]; for RX J1713.7- 
3946 purely leptonic, purely hadronic or also mixed scenarios [7, 12] have been proposed 
until only very recently Fermi-LAT observations firmly shifted the paradigm towards the 
leptonic case. Finally W51C, W28, W49B and IC443 seem to be representative of a class 
of middle-age SNRs associated with MCs, whose TeV emission (not detected in W51C, yet) 
should come from above the cut-off [124]. 

In this respect, observations in the GeV band of TeV-bright SNRs may be extremely 
important for assessing the role of SNRs as hadronic accelerators in that, in this range of ener- 
gies, only relativistic bremsstrahlung and pion decay may produce a photon energy spectrum 
steeper than E~'^. There are in fact a few SNRs, like Tycho and CTB37A, whose GeV-to-TeV 
7-ray emission shows a uniformly steep spectral index: if relativistic bremsstrahlung were 
negligible in the GeV band there would be no lepton-induced mechanisms able to account 
for the broadband emission, therefore strongly supporting a hadronic scenario. In order to 
make such a claim, however, it is necessary to carry out a full non-linear modeling of the 
multi-wavelength emission, from radio to 7-rays, able to self-consistently constrain all the 
different tiles of the observational mosaic. In such an analysis important information may 
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come also by the thermal X-ray emission (both continuous and lines), which can put strong 
constraints on the density of the circumstellar medium and, in turn, on the level of the 
hadronic emission (see [8] for the case of RX J1713. 7-3946, where there are no evidence of 
thermal X-ray emission). 

It is also important to stress that NLDSA calculations may also predict a concave 
spectrum for the accelerated particles, which has to be typically steeper (harder) at the low 
(high) energies: such an evidence of concavity would be of great importance in assessing 
SNRs as very efficient CR accelerators. However, since when the Alfven velocity is taken 
in the amplified magnetic field very large CR efficiencies (larger than 10 — 20%) are neither 
required nor predicted (see also [117]), it is very likely for the deviation from a straight 
power-law from GeV to TeV energies to be well inside any measurement error. 

In any case, since DSA is charge independent, the electron spectrum has to be parallel 
to protons' one, therefore the argument about the necessity of some physical mechanism 
able to account for steep spectra applies also when explaining the 7-ray emission as due to 
relativistic bremsstrahlung, as it has been proposed for Cas A [45]. 

The case of SNR RX J1713. 7-3946 remains however emblematic in highlighting the 
difficulties embedded in the study of a given source: assessing the dominance of one scenario 
over the other required in fact a few years of Fermi-LAT observation. Nevertheless, it is 
important to notice how in figure 3 RX J1713. 7-3946 lies in a quite peculiar position, showing 
a GeV spectrum much flatter than the average: it is likely for RX J1713. 7-3946 not to be 
the most representative 7-ray-bright SNR among the ones listed in table 1. 

It is indeed possible that different emission mechanisms are at work at the same time, 
but in different spatial regions determined, for instance, by a strongly inhomogeneous cir- 
cumstellar medium populated with several MCs. Non-spherical models of such complex 
environments have not been put forward yet, but nevertheless the theoretical framework is 
already developed enough for testing future, detailed, space-resolved maps of the broadband 
emission. In this respect, very interesting information may come in the very next future 
from systems showing emission both from the proper SNR shell and from close MCs, like for 
instance in the cases of W28, W51 or RX J1713. 7-3946. 

5 Caveats and limitations of the present approach 

The simple model put forward in this paper is only meant to provide a first-order quantitative 
insight into the phenomenology of hadronic 7-ray-emission from SNRs. In this section we 
discuss the limitations of the present approach and outline some other physical ingredients 
which may enter a more accurate prediction of 7-ray spectra and fluxes expected from SNRs. 

5.1 Magnetic field modeling 

Both the present treatment of magnetic field amplification and the adoption of an effective 
scattering center velocity are based on a plausible extrapolation in the non-linear regime of 
the quasi-linear theory, and has therefore to be checked against numerical simulations for 
assessing the validity of the assumptions made. We checked a posteriori that the inclusion of 
Bell's non-resonant modes [89, 90], with their standard saturation Pct/Pb — 2c/Kih [91-95], 
does dot change the results above in a sizable way. The very reason of this fact is that Bell's 
mode saturation does not depend on the Alfvenic Mach number, i.e. magnetic field amplifica- 
tion is not expected to become more and more effective the larger the magnetic field is, while 
resonant streaming instability niciy Sctturate to Sj level wliicli eventucilly depends only on 
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(see eq. 2.3). As a consequence, when such an effective saturation is achieved, typical values 
of ^cr make excitation of resonant modes so efficient to overwhelm the possible contribution 
of smaller- wavelength ones. However, the intrinsic complexity and non-linearity of the in- 
terplay between relativistic particles and magnetic fields is such that a deep comprehension 
of how magnetic field amplification really occurs, and hence of how to model it correctly for 
observational purposes, might come only from particle-in-cell simulations (see e.g. [95] and 
references therein). 

It is fair to stress that, from a quantitative point of view, the recipes put forward in 
section 2 have to be regarded as phenomenological tools, in that results might be somewhat 
different if one assumed different prescriptions either for the wave anisotropy, or for the 
effective fraction of the amplified magnetic field entering the Alfven velocity [125] or even for 
the wave transport equation, a solution of which is given by 2.3. 

In the two-fluid model described in section 2 we dealt with relatively small CR accel- 
eration efficiencies, hence the resulting shock are only mildly modified. However, even with 
almost constant velocity and density in the upstream, particles may in principle feel different 
compression ratios because of variations in the local Alfven velocity. A curvature in the 
spectrum may in turn result from the fact that particles with different momenta feel different 
averaged Alfven velocities (i.e. different magnetic fields) in the precursor. Since we assumed 
the magnetic field to be constant upstream as well, no concavity can be consistently pre- 
dicted. The question may be whether this last assumption is realistic or too much simplistic 
and therefore whether spectra might be concave. Strictly speaking, in fact, magnetic field 
amplification should be more efficient where the gradients in the CR distribution function 
are larger, namely in the correspondence of the diffusion length of the particles carrying the 
most of the energy (around 1 GeV for spectra steeper than E""^). Particles with larger ener- 
gies and therefore larger diffusion lengths actually probe regions with lower magnetic field, 
therefore the spectrum may become flatter and flatter at higher and higher energies. 

Some arguments supporting the fact that our two-fluid model can nevertheless be a 
reasonable approximation of a more complex kinetic model are the following: 1) in the far 
precursor non-resonant streaming instability excited by escaping CRs may help to provide 
a substantial magnetic field enhancement, actually reducing the gradient in the magnetic 
field; 2) the possible flattening at the highest energies is however controlled by the largest 
compression factor (r w 4 — 4.5) allowed by ^^ff with ^cr — 0.1, hence the spectral index 
may vary between 2.2-2.5 and 1.9-2 at most. Measuring a possible hint of concavity could 
be very intriguing but present observations, whose systematic + statistical errors on spectral 
slopes are typically in the range 0.1-0.3, allow us only to focus on the more basic question, 
i.e. how it is possible to produce spectra steeper than £'~^ in SNRs. 

Another effect which cannot be easily handled in this simplified approach is the dynam- 
ical role of the amplified magnetic field on the shock dynamics. Pressure in the shape of 
magnetic turbulence may in fact dominate over ordinary thermal plasma pressure upstream, 
leading to non trivial modifications of the Rankine-Hugoniot equations describing the jump 
conditions at the subshock. In the case of resonant Alfven waves such a magnetic feedback 
has been shown to reduce the compressibility of the upstream plasma, providing a smoothen- 
ing of the precursor (i.e. a reduction of the total compression ratio) and, as final consequence, 
steeper spectra of accelerated particles [99, 102, 103]. 

All of these arguments firmly asses the need for the inclusion of a detailed magnetic field 
treatment as a key ingredient in testing hadronic scenarios for the 7-ray emission observed 
from SNRs. Nevertheless, the results outlined in this paper hold as a qualitative description 
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and as a first-order estimates of many of the effects related to tlie observed magnetic field 
amplification ongoing in SNRs. 

5.2 Molecular clouds 

The frequently observed correlation between 7-ray bright regions around SNRs and MCs is 
indeed a strong hint in favor of a hadronic origin for such an emission. There are in fact no 
evident reasons why the most plausible leptonic processes, namely relativistic bremsstrahlung 
and inverse-Compton scattering on the cosmic microwave background or on the Galactic 
infrared radiation, should correlate with the positions of MCs. On the other hand, in a 
hadronic scenario a much more natural correlation is expected because of the huge amount 
of targets provided by environments where the density may easily be hundreds or thousands 
times larger than in the ISM. 

For instance, in order to account for the observed luminosities from objects like IC443, 
W28N, W44, and W51C it is necessary to invoke local density as large as n = 10 — lOOcm"^. 
These sources are associated with MCs and are quite peculiar, in that they show a cut-off 
in the GeV band. They probably represent a somehow distinct class of objects (SNR-MC 
associations?) which cannot be described with the present simple model (they occupy a 
off-set positions in figure 3 as well). Moreover, the luminosity in 7-rays of these middle-age 
sources turns out to be as large as 10^^ — 10'^^ erg s~^, about one or two orders of magnitude 
above the most optimistic predictions in both figure 4 and ref. [2], which typically span the 
range 10^^ — lO^^erg s~^. We cannot but conclude that, at least for these objects, the presence 
of huge repositories of molecular gas is fundamental in order to account for the large fluxes 
observed. 

A detailed analysis of the properties of the 7-ray emission expected from MCs is never- 
theless beyond the goals of this paper, and hence the related phenomenology is not accounted 
for in the formalism above. It is however necessary to spend a few words about their potential 
effects in the matter on debate. 

Regarding the argument about the role of magnetic field amplification put forward in 
section 2, we can make the following considerations. At the zeroth order, the role of a neutral 
component is to boost the luminosity, leaving only the slope unscathed. This is why it may 
be worth focusing more on the spectral slope rather than on the emission level. At the next 
order, one could imagine the shock to slow down significantly because of the larger inertia 
of the swept-up material, but also the magnetic field to be suppressed because of a more 
efficient ion-neutral damping. Since these effects go in opposite directions in changing the 
Alfvenic Mach number, and may partially compensate themselves, it is hard to estimate a 
priori the effects of a relevant neutral component on the particle spectrum. 

The actual role of MCs may however be twofold, depending on whether the forward 
shock has directly crashed into the cloud or not. In the former case, the enhanced density of 
the cloud provides a huge amount of target nucleons, resulting in a magnified, direct probe 
of the CR content of the remnant, while in the latter case a MC may be illuminated only by 
escaping CRs, whose spectrum is not strictly related to the one of particle being accelerated 
at the shock. Some examples of multi-wavelength studies of the 7-ray emission predicted 
from MCs illuminated by CRs escaping from close SNRs are put forward for instance in 
ref. [124, 126-128]. 

In any case, the NLDSA theory needed to predict the spectrum of the particles re- 
sponsible for such an emission has to be completed with an additional ingredient not fully 
understood, yet. In the first case, in fact, one should account for a shock propagating in a 
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partially neutral environment, therefore including the effects of charge-exchange in the hy- 
drodynamics and of ion-neutral damping in the magnetic field treatment; in the second case, 
instead, one should deal with the details of how particle escape from a SNR, a problem related 
with the evolution of the SNR confining power, and in turn, with the evolution of magnetic 
turbulence responsible for particle diffusion. The transition from the diffusive regime in- 
side the SNR and the one in the Galaxy, in fact, spans several order of magnitudes in the 
strength of the diffusion coefficient, D(E); also the energy-dependence of the phenomenon is 
expected to vary from Bohm-like D{E) oc S to D{E) oc E\ with S = 0.3 - 0.6, as inferred 
from secondary to primary ratios in diffuse Galactic CRs. The exact details of such a poorly 
understood transition may dramatically change the predictions for the spectrum of escaping 
CRs impacting into a MC. 

The problem of describing the escape of accelerated particles from a source is particularly 
important also because it is a fundamental piece of information needed to relate particle 
acceleration in SNRs to the diffuse spectrum of Galactic CRs. In ref. [129, 130] the reader 
can find a wider discussion providing some insights into why particle escape from SNRs still 
remains one of the most fundamental open question in the whole theory of the production of 
Galactic CRs. 

5.3 SNR gamma-ray luminosity 

Present estimates of the time evolution of the 7-ray luminosity of SNRs, however, must be 
taken with a grain of salt and always checked against non-linear calculations in which CR 
acceleration efficiency and maximum proton energy are calculated self-consistently. A clear 
example of this is the rather unphysical growth of the predicted emission for very late stages: 
at a certain point the maximum energy of accelerated protons has to drop below a few GeV, 
implying a fading of the hadronic emission. In addition, it is still not clear whether the 
actual release of accelerated particles which have been advected downstream occurs only at 
the "death" of the SNR or whether it is a continuous process, occurring during the whole 
Sedov- Taylor stage [130]. If the latter were the case, it is possible for highest-energy particles 
to leave the remnant well before the beginning of the snowplow stage, hence severely limiting 
the TeV emission in old SNRs. Nevertheless, it is worth remembering that the release of 
high-energy particles from downstream does not affect the global SNR evolution, since in 
our cases most of the energy in CRs responsible for the shell pressure is stored in mildly 
relativistic particles, spectra being steeper than 

The recipe for estimating the total 7-ray emission adopted here and in ref. [2] neglects 
both the time evolution of Emax and the details of the release of advected particles, therefore 
it is usually quite reliable for the GeV emission, but it has to be sometimes regarded only as 
a rule of thumb for multi-TeV photons. 

Another strictly related problem is the calculation of the probability of detecting a 
SNR as a PeVatron, i.e. a source of photons as energetic as a few hundreds TeV produced 
by hadronic interactions of multi-PeV protons. Such an evidence would be, in fact, a clear- 
cut signature for acceleration of hadrons in SNRs, in addition up to the knee observed in 
the diffuse spectrum of Galactic CRs detected at Earth. Again, the estimate of how many 
PeVatrons could be detected cannot fail to self-consistently account for CR acceleration 
efficiency, magnetic field amplification, evolution of Emax and SNR luminosity in 7-rays. All 
of these quantities have to be calculated as a function of the circumstellar medium properties, 
and in particular have to be evaluated at the epoch of the transition between ejecta-dominated 
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and Sedov- Taylor stages. In fact, a SNR might be expected to work as a PeVatron only during 
this peculiar evolutionary stage [113, 130]. 

Another variable that may enter the history of the 7-ray emission from SNRs is the 
abundance of nuclei heavier than Hydrogen, both in the circumstellar medium and in the 
accelerated particles. Taking into account the interstellar abundances of all the chemical 
elements, in fact, increases the rate of nuclear collisions by a factor 1.5-2 [131], while the 
abundances of heavy nuclei in accelerated particles, inferred by the CR chemical composition 
measured at Earth, typically boost the total secondary emission by another factor 2-3, as 
shown in [132]. Very generally, the chemical composition of the wind and of the bubble 
may be quite different from the ISM one, and such a difference might have observational 
consequences on the level of the expected 7-ray emission and, in turn, in the statistical 
considerations above. 

6 Conclusions 

In this work we tried to understand the impact of very recent observations upon the so- 
called hadronic scenario for the 7-ray emission from SNRs. When energetic photons result 
from the decay of neutral pions produced in nuclear interactions between relativistic hadrons 
and background plasma, their spectrum unequivocally maps the one of accelerated particles, 
providing an insight into the instantaneous content of a SNR in terms of cosmic rays. 

In table 1 we collected the most, as far as we know, updated list of 7-ray-bright SNRs 
and of 7-ray-sources potentially associated with SNRs, showing that basically all of their 
energy spectra are invariably steeper than E~'^. As outlined in section 1, such steep spectra 
are apparently at odds with standard expectations of NLDSA theories at SNR shocks, which 
predict the spectra of high-energy particles to be flatter than E~'^ when the acceleration is 
efficient. 

In section 2 we illustrated a possible way to reconcile NLDSA predictions with observa- 
tions, namely by accounting for the generation of magnetic turbulence induced by the same 
accelerated particles. In fact, when resonant streaming instability is as effective as to re- 
produce the large magnetic field inferred in young SNRs (a few hundreds /iG), the standard 
picture of the particle-wave scattering is, very likely, no longer suitable for describing CR 
transport. In particular, the velocity of the magnetic perturbations particles scatter against 
may be not negligible with respect to the fluid velocity so that, when this occurs, a sizable 
modification of the average compression ratios felt by diffusing particles is expected, even- 
tually affecting the spectrum of accelerated particles. Such an effect has already been put 
forward before (see e.g. [19]), but its quantitative implications, especially under the light of 
the more recent discovery of effective magnetic field amplification in young SNRs, have not 
been deeply investigated yet. This key ingredient may actually reverse the role of the CR 
acceleration efficiency in the determination of the spectral slope, predicting that the larger 
the CR pressure, the more efficient the amplification of the magnetic field and, in turn, the 
steeper the spectra of the accelerated particles. 

It is also worth stressing that the need for a physical mechanism producing steep spectra 
is required when explaining the detected emission as either due to relativistic bremsstrahlung 
or to pion decay, since in both cases the produced 7-ray spectrum is parallel to the energy 
distribution of the parent particles. Strictly speaking, only inverse-Compton emission would 
not require to accelerate particles with a spectrum steeper than E~'^ but, on the other 
hand, such an emission mechanism tends to produce rather flat photon spectra, at odds with 
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broadband observations listed in table 1 (see also comments in section 4.2). The spirit of the 
calculations put forward in section 2 may thus apply not only in hadronic scenarios, but even 
in several other situations where the emission is thought be predominantly leptonic and due 
to relativistic bremsstrahlung, like in Cas A. 

We showed in section 4 that a CR acceleration efficiency between 10 and 20% may 
easily account for the spectra observed in most of the detected SNRs (figure 3). In our 
calculations, the expected 7-ray spectrum is calculated during the different evolutionary 
stages of a remnant produced by a core-collapse SN. This means that we also included the 
non trivial circumstellar profile induced by stellar winds launched during the red-giant and 
Wolf-Rayet stages of the SN progenitor. In addition, we worked out the SNR luminosity in 
the GeV (namely, above 100 MeV) and in the TeV band and found that it has a minimum 
when the remnant is between 5 and 10 thousands years old, corresponding to the propagation 
of the forward shock into the hot and rarefied bubble excavated by the Wolf-Rayet wind, a 
feature that seems to be recovered also in figure 3. 

Finally, in section 5 we discussed strong points and limitations of the present approach, 
highlighting open theoretical problems, as outlining a consistent description of the non-linear 
interplay between CRs and magnetic field, but also observational issues, as the observed 
correlation between strong 7-ray emissivity and MCs. 

Our findings, worked out in the context of a rather phenomenological scenario, have 
however to be checked against fully non-linear calculations able to self-consistently account 
for shock dynamics, particle acceleration and magnetic field amplification, but they indeed 
represent a preliminary, necessary step for understanding hadronic emission from SNRs. 
Under the light of this revised theoretical framework for the hadronic emission, we finally 
discussed what future observations achievable with the present (and next) generation of 7- 
ray telescopes may tell us, not only regarding non-thermal properties of SNRs, but also 
regarding the origin of Galactic CRs. 
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7-ray-bright Supernova Remnants 
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[14, 34, 44] 
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7 Cygni 
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~2.17 
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" W28N, IC443, W51C, W44 and W49B show evidence of a cut-off around 1-20 GeV. 
Ref. [14] points out SNRs in physical contact with MCs, except Tycho according to ref. [78]. 



Table 1. Photon spectral index a inferred in 7-ray-bright SNRs, both in the GeV and in 
the TeV bands. Associated systematic errors are typically in the range ±(0.1-0.2). Informa- 
tion about SNR nomenclature, ages, distances and associations can be found in the Green's Cat- 
alogue (http://www.mrao.cam.ac.uk/surveys/snrs/snrs.data.html) and in refs. [79-81]. See also 
http://tevcat.uchicago.edu and http://www.mppmu.mpg.de/~rwagner/sources/index.html for cat- 
alogues of known TeV sources. 
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